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Введение. Принято считать, что электроупругие модули 
не зависят от амплитуды и частоты колебаний. Это под- 
ход отражен в российских и зарубежных стандартах для 
определения полного набора электроупругих модулей 
пьезокерамики. Например, для определения пьезомодуля 
Яз1 образца в форме диска необходимо провести изме- 
рения в трех частотных областях: в области первого и 
второго резонансов, в области антирезонанса и на часто- 
тах много ниже 1 кГц. В связи с этим предполагается, 
что при определении 431 модули исследуемой керамики 
в диапазоне частот от 1 КГц до второго резонанса не 
зависят от частоты. 

Целью данной работы является исследование зависимо- 
сти от частоты электроупругих модулей керамики. При 
этом используется образец в форме диска из ЦТС (цир- 
коната-титаната свинца). 

Материалы и методы. Использованы методы постанов- 
ки и решения задач стационарной электроупругости и 
разделы теоретических основ электротехники. Для реа- 
лизации метода конечных элементов применены метод 
возмущений и пакет прикладных программ АМЗУ5. Ре- 
зультаты экспериментов обработаны в среде МАТГАВ. 
Результаты исследования. Для пьезокерамики ЦТС 
исследованы зависимости от частоты различных моду- 
лей (пьезоэлектрических 431, диэлектрических =1. и 
упругих модулей гибкости $11, 512, 513). Рассмотрены 
радиальные колебания образца в форме диска с электро- 
дами на торцах. Толщина образца — 1 мм, диаметр — 
40 мм, диапазон колебаний — до 700 КГц. Сначала ис- 
следовалась частотная зависимость для упругих модулей 
керамики из определения десяти резонансных частот. 
Затем по измеренным значениям проводимости образца 
была определена зависимость от частоты модулей 
431 И #13. С этой целью использовалось полученное в 
работе выражение для электрической проводимости из 
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решения радиальных колебаний диска с учетом его тол- 
щины. 

Обсуждение и заключения. Разработана методика для 
определения зависимости от частоты модулей пьезоэлек- 
трической керамики ЦТС. Образец в форме диска иссле- 
довался при диапазоне частот 15-650 КГц. Показано, что 
в диапазоне до 650 КГц упругие модули $11, $12, 513 © 
верхним индексом ЕЁ (он опущен) или измеренные при 
постоянном электрическом поле практически не зависят 
от частоты. В указанном диапазоне константы #23, 41, ре 
для рассматриваемых радиальных колебаний имеют не- 
значительную частотную зависимость. 


Ключевые слова: пьезоэлектрическая керамика, элект- 
роупругие модули, схема замещения, зависимость от 
частоты, диск, метод конечных элементов, АМЗУ5, ме- 
тод возмущений. 
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Введение. Значительное число работ посвящено исследованию и развитию математических методов 
решения задач о колебаниях пьезоэлектрических тел. В составе пьезоэлектрических устройств пьезоэлементы 
служат для возбуждения и регистрации колебаний, вызванных внешними воздействиями. Выбор пьезоэлектри- 
ческого материала для измерительных преобразователей и анализ их характеристик требует большого объема 
сведений о параметрах материалов. К таким сведениям относятся: 

— полный набор электроупругих модулей [1], 
— потери, 
— эквивалентные электрические схемы, или схемы замещения пьезоэлементов [2]. 

Реактивные динамические параметры Г, С схем замещения определяются упругими диэлектрическими 
и пьезоэлектрическими модулями, а также плотностью пьезокерамики [3]. Принято считать, что электроупру- 
гие модули — параметры схем замещения постоянны, они не зависят от амплитуды (слабые электрические по- 
ля) и от частоты колебаний. Все это нашло отражение в российских [4] и зарубежных [5, 6] стандартах для 
определения полного набора электроупругих модулей пьезокерамики. Например, для определения пьезомодуля 
Яз1 образца в форме диска необходимо провести измерения в трех частотных областях: в области первого и 
второго резонансов, в области антирезонанса и на частотах много ниже резонансов 1 кГц. В связи с этим пред- 
полагается, что при определении 431 модули исследуемой керамики не зависят от частоты в достаточно широ- 
ком диапазоне от 1 КГц до второго резонанса. 

Целью данной работы является исследование зависимости от частоты электроупругих модулей керами- 
ки. Рассматривается образец в форме диска из ЦТС (цирконата-титаната свинца), т. к. именно такой тип кера- 
мики наиболее хорошо известен. Для проверки корректности представленной методики используется програм- 
ма АМЗУ$ [7]. 

Материалы и методы. Рассмотрим пьезоэлектрический диск толщиной 2 и радиусом а. Введем ци- 
линдрическую систему координат (и, ©, 2), причем ось 2 совпадает с направлением оси поляризации. Коорди- 
натная плоскость 2 = 0 совпадает со срединной плоскостью диска. 

Основываясь на известных линейных пьезоэлектрических соотношениях, уравнениях динамики 
сплошной среды [1] и Максвелла [8], систему уравнений осесимметричных колебаний пьезоэлектрического 
диска можно записать в следующем виде: 


И + роз? = 0, 


.Т,„ + 051, + + ро? =0; д.0, +050, += =0. (р 
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Здесь и далее введены следующие обозначения и определения: И, Й’ — механические смещения по координат- 
НЫМ ОСЯМ Г, 2 соответственно; © — круговая частота; р — плотность; Ти„ — механические напряжения; 2.., О» 
— компоненты вектора электрической индукции; д/ и дз — операторы производных по ги 2. 

В случае осевой поляризации линейные уравнения пьзоэффекта для слабых электрических полей в ци- 
линдрических координатах можно записать [1]: 


Тек = с11 010 + с12 =  с1здз И + ез1д3 Ф, 

Тев = с1204 И + сл - + с1з03 И/ + ез1 03 Ф, 

Т.» = с1з910 + С1з- ой Сзздз И/ + езздз Ф, (2) 
р, = ез1 (010 + =)  езздз И/ — =. ф; О, = е15(0зИ + д.И/) — Е11 д: Ф. 


В соотношениях (2) и далее введены следующие обозначения: сии — упругие постоянные в матричном 
обозначении, измеренные на образцах с закороченными электродами или при постоянном электрическом поле 
Е (верхний индекс Е у констант керамики здесь и далее опущен); еи„ — пьезоконстанты; =и„ — диэлектриче- 
ские проницаемости для зажатого образца; ф — электрический потенциал, причем Ё = —ота4 ф [2, 8]. 

Будем считать, что на электродированных торцах и на боковых поверхностях диска заданы граничные 
условия [1, 2]: 











при 2 ==й Т,„„ = 0; Т,, = 0; ф=2И, (3) 
приг=а Т,., = 0; Т,.. =0; О, =0. (4) 
В (3) величина 27 — подведенная к торцам электрическая разность потенциалов [2, 8]. 
Введем безразмерные координаты и величины по формулам: 


г 2 в р 13 
=-; ЕЕ; ==; О = ов |2; 0. = ха а= 11 = —^ 
5 ые в’ а: Сад’ 11 сз 


Решение краевой задачи (1-3) состоит из суммы двух решений: 
— однородного решения при нулевых граничных условиях при 2 = -Й; 











— частного решения, удовлетворяющего только ненулевым условиям на торцах (3). 

Построенная система однородных решений позволит в дальнейшем удовлетворить и граничным усло- 
виям (4) на боковой поверхности (как правило, с помощью вариационных методов). 

Не составляет труда построить частное решение 02 = соп5{ и 08 = 0, которое автоматически удовле- 
творяет третьему уравнению из системы (1) и граничным условиям на торцах (3). Механические и электриче- 
ские составляющие частного решения согласно первым двум уравнениям системы (1) равны: 


ЗС 
То = ТФ = Або +) изтОО + ен; И° = 0; 1% =0; 
33 
фо = КК + АЙзт(ХО; 02 = —ЕззК; То = еззК + ез.Ах соз(хО; р 
у р у 1 © 


е вы | 
В езз (0500 — КЁ 51100)’ 1 _ 2 ав 00 
Хх 





В (5) и далее введены следующие обозначения: 


Е С33Е с 
с3з 


с3з' 





Если вектор внешних усилий и о потенциал равны нулю на торцевых поверхностях, то 
построение однородных решений связано с определением корней дисперсионного уравнения [9]. Для симмет- 
ричных колебаний дисперсионное уравнение имеет вид: 

а„М, вап 1(Ви) = 0, (и=1,2, 3). (6) 

В (6) введены следующие обозначения: 

ав = а? с1зК1п + СззК2п + @ззКзп; Ви = Киви + Кап; 
М, = Б-Кзз — ВзКз2; М» = БзКз1 — В1Кзз; Мз = В1Кз2 — Б>ЁКза. 
Здесь Кии — алгебраические дополнения элементов третьей строки определителя системы для симметричных 
колебаний (1); В„ — корень бикубического уравнения из [9], и — безразмерное волновое число колебаний по 
ОСИ 7. 
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Для нахождения корней о при заданных значениях © необходимо совместно решить дисперсионное (6) 
и бикубическое уравнения. Подробный анализ корней дисперсионного уравнения симметричных колебаний без 
потерь для пьезоэлектрического слоя приведен в [9], с учетом потерь — в [10]. Именно сумма частного и одно- 
родного решений позволит удовлетворить граничным условиям как на торцах, так и на боковых поверхностях 
диска. 

При произвольном соотношении размеров диска обратная задача о его вынужденных колебаниях (1-4) 
является очень сложной для анализа и имеет конечное аналитическое решение лишь в некоторых частных слу- 
чаях (например, колебания тонкого диска или без учета его толщины, когда # «< 1). Поэтому при определении 
модулей керамики обратную краевую задачу удобнее решать приближенными методами с учетом толщинных 
поправок для низкочастотных радиальных колебаний диска. В этом случае можно получить аналитическое ре- 
шение в виде конечных формул. В данной работе ищем решение в виде разложения по малому параметру Е: 

42 = =202 + 7=* 0% + 1=8 06... Ц =Ц + #2020), + =040, -.... (7) 








Здесь И„ — вектор-функция с составляющими И(0, И’, ф); п — порядок построенной приближенной теории 
рассматриваемых в работе симметричных колебаний диска; 7, | — неизвестные постоянные, которые зависят 
от модулей пьезокерамики, определяются из (1, 2) и удовлетворяют однородным (нулевым) граничным услови- 
ям (3). 

Опускаем громоздкие выкладки. Ограничимся членами с =6 для определения волнового числа а, чле- 
нами с =? для вектор-функции И в (7) и приведем конечный результат рассматриваемой краевой задачи с вы- 
полненными граничными условиями (3) на торцах [11]: 











у= С: , — 912. 10 2431 (в _ 613 _ `` 913. 
Зе: ° 51’ Г 23 (5аа + 512), Сз3 511 +512, 
о г +у] +7: 
=у|— 1е20; 
ПУ 15+ 5. У] +Р 
12 аз1 





а о | 
{е20 = — с. (24:1 +а . 
р 45СР. Е13(1-К2) 14912 11( 31 з3)| 


Здесь 5„„ — модули гибкости при постоянном Е; 431 — пьезомодуль; К» — планарный коэффициент связи; у — 
планарный коэффициент Пуассона; =13 — свободная диэлектрическая проницаемость диска. 

Далее введем следующие определения и обозначения: &5 = =202 — приближенное волновое число ну- 
левого порядка; 42 =Е?202 +1='0% — —  приближенное волновое число второго — порядка; 
02 = =20 7+ =“ + 1=б Од ° — приближенное волновое число четвертого порядка; С — емкость диска. 

В табл. 1 представлены результаты точного решения волнового числа о керамики РИТА в зависимости 








от частоты из дисперсионного уравнения (6) и приближенный расчет из (7) при # = 0,033 для различных при- 
ближений о. 

















Пример расчета волнового числа а при различных частотах для рассматриваемого а —— 
КГц о а, из (6) 60 02 94 с р1его = 0 04 
50 0,085 0,045055 0,045049 0,045055 0,045055 0,045055 
250 0,425 0,2262 0,2252 0,22603 0,22605 0,2261 
500 0,85 0,4578 0,4505 0,4568 0,4574 0,4575 
700 1,1903 0,6522 0,6307 0,64777 0,6514 0,6518 





























Приведенные в табл. | результаты показывают, что для диска толщиной | мм и частот до 700 КГц 
расчет волновых чисел из дисперсионного уравнения (6) и приближенный расчет для 04 практически совпада- 
ют, а поправкой рЁего можно пренебречь. В этом случае разложение (7) для а? практически зависит только от 
геометрических размеров диска, плотности и модулей гибкости $11, 512, 513. 

Опуская относительно громоздкие выкладки, приведем выражение проводимости У для пьезоэлектри- 


ческого диска. Оно получено из приближенного решения (7) для различных =. 
7000 — нулевое приближение # = 0, или известное уравнение радиальных колебаний диска нулевой толщины: 


зы 96: оббы Воее 
Е. В и хо (х) — 1-1 (х)/ 


7040 — приближенное волновое число 4-го порядка: 
ча. _ И. ло 
х=-; У040 = «С (1 Кр + Кр (1 + у) ны) (8) 


7042 — приближенное волновое число четвертого порядка, частное решение второго порядка: 


Мадорский В. В. и др. Исследование зависимости от частоты констант поляризованной пьезокерамики в схемах замещения 


МааогзКу И. У. 1 а. 5шау оп гедиепсу 4ереп4епсе оГроатсе4 рёезосегат!с$ сопуатз т едшущет сдгсий5 








02е2. х 
нова = с (1 Вена Л ) 


З= 333 Сзз О Ло(х) — (1 — у)/1 (<) 


7242 — приближенное волновое число четвертого порядка, однородное и частное решение второго порядка: 





02 е2. И ) т @ззС1з 
узна = ес (1 Чена + Ка + У) (1 тих 2 ›) 
2зз с3з е31С3з 


1: 
эт = Боб — а - Уре [1 + 208 (5+ : = аду хо 69; (9) 


ЕззС1з _ 431 
езз 511 + 512 
В табл. 2 представлены результаты численно-аналитических расчетов. Проводимости точного реше- 





Р _ 
ез1 = ез. + 


ния (5, 6) и приближенных решений (8, 9) по =2 найдены с помощью программы МАТЕАВ [12]. Для диска из 
РИТА [13] толщиной | мм и диаметром 30мм численные расчеты проводились в системах АМЗУЗ и 
АСЕГАМ [14]. Последний пакет программ разработан сотрудниками Южного федерального университета 
(ЮФУ) и ориентирован на расчеты пьезоэлектрических устройств. 














Таблица 2 
Расчет проводимостей для различных частот 
Частота, КГц |АМЗУ5, См/м |Решение (5,6), См/м | 7000, См/м | 7040, См/м |042, См/м | 7242, См/м 
350 5,03 5,032 4,9688 5,014 5,0457 5,028 
554,4 8,74169 8,741029 8,478 8,5864 8,73237 8,735 
700 10,98 10,97 10,376 10,682 10,957 10,96 





























Из табл. 2 видно, что при построении приближенных теорий типа (8, 9) желательно использовать раз- 
ложение (7) для расчета волнового числа как минимум второго порядка. 

Результаты исследования 

1. Исследование зависимости упругих модулей керамики от частоты. Для достаточно хорошо изу- 
ченной пьезокерамики ЦТС проведем исследование зависимости от частоты следующих модулей: пьезоэлек- 
трических 431, диэлектрических #13 и модулей гибкости $11, $12, 51з. Рассмотрим низкочастотные радиальные 
колебания образца в форме диска с электродами на торцах. Толщина образца — 1 мм, диаметр — 40 мм, диапа- 
зон колебаний — до 700 КГц. Сначала исследуем частотную зависимость для упругих модулей. Для этого, со- 
гласно [15], предварительно измеряем три первые резонансные низкие частоты {. (основной резонанс и два его 
обертона). Упругие константы определяют из решения частотного уравнения (9) радиальных колебаний диска с 
учетом относительной толщины : три уравнения для трех неизвестных Именно относительная толщина # отли- 
чает приведенное выше частотное уравнение (9) для радиальных колебаний диска конечной толщины от из- 
вестного частотного уравнения для радиальных колебаний диска с «нулевой толщиной» [15]: 

В.Л(®) = 1-У) л(®. 

Введение толщинной поправки & в решение известного уравнения радиальных колебаний образца в ви- 
де диска повышает точность и информативность измерения упругих модулей. 

Определенные по [11, 15] методом трех резонансов упругие податливости для рассматриваемой кера- 
мики оказались равны: 

511 = 12,29е- 12, 512 = -4,05е- 12, 53 =—5,28е - 12. 
В табл. 3 приведены первые десять резонансных частот для керамики ЦТС19. В стоке «Эксперимент» 





указаны частоты, измеренные на установке для измерения проводимости И’К 6500В в Институте высоких тех- 
нологий и пьезотехники при ЮФУ. В строке «Расчет» — рассчитанные по формуле (9) для независимых от ча- 
стоты упругих модулей $511, 512, 51з, определенных по [11, 15]. Погрешности измерений не превышали значений, 
рекомендуемых стандартом [4]. 











Таблица 3 
Первые десять резонансных частот для керамики ЦТС19 
Резонансные частоты /,, КГц 
Эксперимент 50,8 132,8 |211 2884 |364,56 |4394 |513,2 |5854 |6544 |722,8 
Расчет 50,8 132,8 |210,95 |288,1 |364,8 |440 |516 590,4 |663 735,4 









































Из табл. 3 видно, что в диапазоне частот до 600 КГц различие между расчетными и экспериментальны- 
ми данными при постоянных упругих модулях не превышает 1 %, поэтому можно сделать следующий вывод: 
«Для рассматриваемой керамики упругие модули $11, $12, $1з с верхним индексом Ё (он опущен) или измерен- 
ные при постоянном электрическом поле практически не зависят от частоты в диапазоне до 650 КГц». 
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2. Исследование зависимости от частоты модулей 431 и #213. Используем выражение (8) для прово- 
димости У040 из решения радиальных колебаний диска с учетом его толщины, а также измеренные при ком- 
натной температуре значения мнимых частей комплексной проводимости исследуемой пьезокерамики. В этом 
случае можно исследовать зависимость модулей 431 и =13 от частоты в диапазоне от 10 КГц до 600 КГц. Для 
определения двух неизвестных 41 и #1. измерим У при двух значениях частот /, №. Разница между № и Л вы- 
брана так, что можно пренебречь зависимостью модулей 41 и =1з от частоты в диапазоне /- р. 

В настоящей работе введены два условия: 

 Б-Лл=200 Гц; 

2) проводимости рассчитаны для идеального пьезоэлектрического тела или без учета потерь, поэтому 
не учитываются частоты резонансов или в окрестности резонансов. 

В итоге получим легко решаемую систему двух линейных уравнений относительно двух неизвестных 

2: Г 
Кри 233. 

На рис. 1 представлены зависимости К ‚Чи 213 /Е0 от частоты 10-650 КГц. 
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Рис. 1. Зависимости пьезоэлектрических параметров от частоты 


Здесь & — диэлектрическая проницаемость вакуума. В диапазоне от 15 до 650 КГц коэффициент связи 
| сначала увеличивается с ростом частоты от 0,28 до 0,34, а затем убывает до 0,26. Пьезомодуль 431 имеет ана- 
логичную частотную зависимость: при 15 КГц с —1,50е '° до максимума 1,661е '°, а затем убывает до 1,31е ®. Отно- 
сительная диэлектрическая проницаемость с ростом частоты монотонно уменьшается с 1766 до 1455. 

На рис. 2 представлены зависимости от частоты в диапазоне 1-10 КГц измеренной и рассчитанной по 
формуле (9) емкости керамики при постоянных модулях, определенных согласно ОСТ [4, 15]. 
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Рис. 2. Зависимость емкости пьезокерамики от частоты 
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Видно резкое отличие двух емкостей керамики ЦТС для низких частот. Вероятно, этим и объясняется 
зависимость параметров керамики от частоты ниже 15 КГц (см. рис. 1). Это предполагает наличие более высо- 
кой степени дисперсности — в первую очередь, у пьезомодуля 431. Подробнее этот вопрос будет рассмотрен в 
следующей статье (константы #зз^Т, 41, А» исследуются в низкочастотном диапазоне). 

Обсуждение и заключения. Разработана методика для определения зависимости от частоты модулей 
пьезоэлектрической керамики ЦТС. Образец в форме диска исследовался при диапазоне частот 15-700 КГц. 
Показано, что в диапазоне до 650 КГц упругие модули $11, 512, 51з с верхним индексом Е (он опущен) или изме- 
ренные при постоянном электрическом поле практически не зависят от частоты. В указанном диапазоне кон- 
станты 233, 41, К» для рассматриваемых радиальных колебаний имеют незначительную частотную зависимость. 

Изучена зависимость только реальной части модулей керамики от частоты, потери не учитывались. По- 
этому экспериментально полученные зависимости мнимой части проводимостей керамики от частоты были 
измерены вдали от резонансов, где влияние потерь ничтожно. Тем не менее в дальнейшем необходимо учиты- 
вать как реальные, так и мнимые части модулей. Это значит, что потери являются включенными. Данную про- 
блему предполагается изучить с учетом зависимости комплексных модулей исследуемой керамики от частоты. 

В представленной работе рассмотрены радиальные колебания образца для получения более полной ин- 
формации при измерении большего набора констант (кроме пьезоэлектрических и диэлектрических модулей 
измерены упругие модули $11, 512, 51з ). Данный метод может быть распространен и на другие формы образцов 
(стержни, пластины и т. п.), поскольку в этих одномерных задачах происходит простая замена функций Бесселя 
на тригонометрические функции. 
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